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Abstract 

In the context of a semiclassical approach where vectorial gauge fields can be considered 
as classical fields, we obtain exact static solutions of the SU(N) Yang-Mills equations in a 
(to + 1) dimensional curved space-time, for the cases to = 1,2,3. As an application of the 
results obtained for the case to = 3, we consider the solutions for the anti-de Sitter and 
Schwarzschild metrics. We show that these solutions have a confining behavior and can be 
considered as a first step in the study of the corrections of the spectra of quarkonia in a 
curved background. Since the solutions that we find in this work are valid also for the group 
£7(1), the case to = 2 is a description of the (2 + 1) electrodynamics in presence of a point 
charge. For this case, the solution has a confining behavior and can be considered as an 
application of the planar electrodynamics in a curved space-time. Finally we find that the 
solution for the case to = 1 is invariant under a parity transformation and has the form of a 
linear confining solution. 

Keywords: Yang-Mills equations, exact static solutions, semiclassical approach, curved 
space-time. 



1. Introduction 

The quark confinement problem [1, 2, 3, 4, 5] cannot be solved by the perturbation 
theory of quantum chromodynamics (QCD) so long as the confinement is a nonperturbative 
phenomenon [6]. In essence the question is about how to confine the colour charges at 
distances of order of characteristic hadronic sizes [7]. The semiclassical approach permits 
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to acquire non-perturbative information about QCD starting from the solutions of classical 
partial differential equations of SU(3) Yang-Mills theory [8, 9]. The searching of solutions 
of classical Yang-Mills equations in presence of static external sources has been motivated 
by understanding the role of non-abelian gauge fields of quark confinement problem. One 
of the first works about this subject showed that if an external source is distributed over a 
thin spherical shell the Coulomb solution is unstable in a specific regime [7]. 

A wide range of solutions of (3 + l)-dimensional Yang-Mills equations in presence of 
localized and extended external sources can be found in the literature [8]- [33]. Likewise some 
exact retarded solutions to Yang-Mills equations with sources composed of iV arbitrarily 
moving coloured point particles were studied in [34]. On the other hand some specific 
solutions of the (2+l)-dimensional SU (2) Yang-Mills equations were obtained in [35, 36]. The 
discovery of global regular solutions of SU (2) Einstein- Yang-Mills equations [37] originated a 
great interest about spherical symmetric solutions [38]- [43]. For instance, the solutions of the 
reduced SU(2) Einstein- Yang-Mills equations with spherical symmetry was presented in [44]. 
Additionally some spherical solutions were considered in [45] for the SU(2) Einstein- Yang- 
Mills theory with a negative cosmo logical constant. Concurrently, it was given a rigorous 
proof for the existence of infinitely many black hole solutions to the 577(3) Einstein- Yang- 
Mills equations [46]. Latterly, the study of quark confinement problem in curved space-time 
has been a subject of interest [47]. 

A semiclassical approach motivated by the black hole physics techniques was proposed 
by Goncharov [48, 49, 50] to describe the energy spectra of quarkonia by solving the Dirac 
equation in presence of SU(3) Yang-Mills fields representing gluonic fields. In the context 
of this approach, explicit calculations have shown that gluon concentration is huge at scales 
of the order of 1 fm [51]. The solutions obtained there can model the quark confinement 
satisfactorily, suggesting that its mechanism might occur within the framework of QCD [52]. 
This implies that at large distances the gluons form a boson condensate and, therefore, gluons 
can be described by the classical S"[/(3)-gauge fields [52]. For this reason the dynamics of 
the strong interaction at long distances could be described by equations of motion of the 
SU(3) Yang-Mills theory [52]. 

In this work we follow the semiclassical approach proposed by Goncharov, presented in 
detail in [52], and we obtain some exact static solutions of the (n + l)-dimensional SU(N) 
Yang-Mills equations in curved space-time, for the cases n — 1,2, 3. As an application of 
the results obtained for the case n = 3, we consider the solutions for the anti-de Sitter and 
Schwarzschild metrics. We show that these solutions have a confining behavior and can be 
considered as a first step in the study of the corrections of the spectra of quarkonia in a 
curved background. Since the solutions that we find in this work are valid also for the group 
£7(1), the case n = 2 is a description of the 2 + 1 electrodynamics in presence of a point 
charge. For this case, the solution has a confining behavior and can be considered as an 
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application of the planar electrodynamics in curved space-time. i.e. electrodynamics in two 
spatial dimensions (QED2+1). For the case n = 1, we find that the solution is invariant 
under a parity transformation and have the form of a linear confining solutions. 

The structure of this paper is as follows. In section 2 we first present some preliminary 
aspects related to basic definitions and notation. In section 3, we obtain the solutions for the 
(n+l)-dimensional SU (3) Yang-Mills equations in curved space-time for the cases n = 1,2, 3. 
Finally our conclusions are summarized in section 4. 



2. Preliminaries 

We work on some spacetime manifold M where in local coordinates the line element looks 

as 

ds 2 = g fiu dx tl dx u , (1) 

and the components g^ v take different values depending on the choice of coordinates and 
dimensions. The Hodge start operator * is defined by relation: A P (M) — > A n ~ p (M), where 
A P (M) is the space of p-form over the manifold M under consideration. If {dx 1 , dx n } is 
the base for A P (M) then 

(l/2) 

*(cte ?1 A ... A dx 1 *) = * _ g' 1 ' 1 ...^^,,,^,,,^^ 1 A ... A dx 1 ". (2) 

The exterior differential, which is represented by d, is defined as: A P (M) — > A P+1 (M). This 
means that 

d = d„dx". (3) 
If the connection A in the gauge group SU(N) is defined as 

A = A^dx" = A^XJx", (4) 

where A a ^ are the non-abelian fields associated to SU(N), \ a are the generators of SU(N), 
with a = 1, 2, N 2 — 1, then the curvature F can be defined using the exterior differential 

as 

F = dA + gAAA = F^Xadx 11 Adx\ (5) 

where F^ u represents the non-abelian stress tensor associated to the gauge group SU(N). 

It is possible to write the SU(N) Yang-Mills equations of the non-abelian gauge field in 
presence of sources using the Hodge star operator as follows 

d* F = g(*F AA- AA*F) + gj, (6) 
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where F is the curvature (5), A is the connection (4), g is the gauge coupling constant 
associated with the gauge group SU(N) and, for example, if considering the QCD-lagrangian 
corresponding to that group then J is a nonabelian current given by 

J = j;\ a * (dx») = *j = *{]«X a dx») = $(/ ® lll )\ a ^\ a dx^ (7) 

where \& are the Dirac fields, 7° are matrices which represent the Clifford algebra and I 
represents the identity. For the case of a point particle at rest, this current density is given 
by J = *(j^X a dx^) = 5(r)q a X a * dt, where q a are constants and then q a X a = T is a constant. 
In what follows we put J = 5{r)q a X a * dt. 

In an analogous way as it is performed into the functional quantization procedure of 
Yang-Mills theories, we fix the gauge through the condition 

div(A) = -U M (V&TA0 = 0. (8) 

3. Solutions in a (n + l)-dimensional curved space-time 

The SU(N) Yang-Mills equations are a non-linear system of coupled partial differential 
equations. In this section we will present some static and exact solutions of the SU(N) 
Yang-Mills equations in a (n + 1) dimensional curved space-time, for the cases n — 1, 2, 3. 

3.1. Case n = 3 

The metric for a curved static space-time and spherical symmetry is specified by 

ds 2 = g^dx^dx" = a 2 {r)dt 2 - (3 2 {r)dr 2 - r 2 [dQ 2 + sin 2 6d^ 2 ). (9) 

Due to we are interesting in the description of problems which involve spherical symmetry 
in which the charge point is localized in r = 0, then we assume that the connection (4) has 
the following functional dependence 

A = A^{r)dx^ = A t (r)dt + A r (r)dr + A e {r)d6 + A v (r)d(p. (10) 

For this case the gauge condition (8) allow to the following equation 

d r ( r2a( /3 r ( ) r ^ (r) ) + a(r)(3(r)A e (r)cot9 = 0, (11) 

that can be satisfied if Ag(r) = and A r (r) = ■ We can set C = in the A r solution 

because this fact does not affect the form of A t and A^ solutions. So we assume the following 
anzats for the form of these solutions [53] 

At =/(r)r, (12) 

A v =g(r)A, (13) 
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where T and A are linear combinations of the group generators. As the exterior differential 
in spherical coordinates is 

d = d t dt + d r dr + d e dQ + d^dip, (14) 

then the curvature is given by 

F = dA + gA A A = -d r f(r)Tdt A dr + d r g(r)Adr A dip + gf{r)g(r)\T, A]dt A dip. (15) 

For this case the SU(N) Yang-Mills are written by the equations (6), but with the connection 
given by (10) and the curvature by (15). Applying the Hodge start operator to F with 
inserting the result into (6) and using the following relations 

StTlO 

*(dtAdr) = d9 A dip, 

a(r)p(r) 

*(dt A d9) = ^^\ Sm ° dr A dip, 
a{r) 

-0(r) 

*(dt A dip) = -. \ K ' dr A d6, 
a(r)sznO 

— ot( r) 

*(dr A dip) = - , ' dt A dd, 
p(r)szn9 

a(r)sin9 , 
*(dr A dff) = ' , dt A dip, 
P{r) 

*(deAd<p) = °%^- dthdr, 
r 2 sin0 

we can obtain, for r ^ 0, that the SU (N) Yang-Mills equations are reduced to the following 
coupled equation system 

dr ^^y = - m £- )f (rr 9 (r)[r,[r,Al (16) 

^( ^W e^T = 9 2 f(r)9(rmr) + rV(r)g(r) y ^ Q 

\a(r)p{r) J a(r) a(r) 

A non-trivial solution from this equation can be obtained if the coupled equations (16) 
and (17) satisfy the abelian condition [A, [r, A]] = 0, then we obtain the two following 
independent equations 

( a(r)d r g(r) \ 

H W) ) = ' (18) 

( r 2 d r f(r) \ r 2 /3(0)^(r) 
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that have the following solutions 



g(r) = 

f(r) = 

We note that the abelian condition [r, A] = is satisfied in a non-trivial way if and only if 
one of the following conditions hold: 

1) If within the combinations of T and A, in terms of the generators of the group, it is 
present only the matrix that constitutes the Cartan subalgebra of the SU(N)-Lie algebra 
i.e. a maximal abelian subalgebra. This means that the commutator of any two matrices of 
Cartan subalgebra is equal to zero. 

2) If T = kA and k is a constant. 

C.l. Case of flat space-time 

We can obtain from (20) and (21) that the solutions, for the flat space-time case a = 
= 1, are 

g(r) = b 1 r + B 1 , f(r) = -a 1 /r + A 1 . (22) 

These solutions are confining solutions in a semiclassical approach where the gluons can 
be considered as classical fields [48, 49, 50, 51, 52, 54, 55, 56, 57]. We have obtained the 
confining solutions (22) using the ansatz given by (12) and (13). We note that these solutions 
were obtained in [49], but using the ansatz 

A = r^a a X a dt + A r dr + r»"f\ a d6 + r Va a X a dip, (23) 

where \i a , u a , p a , a a , ■y a and (3 a are arbitrary real constants. The confining potentials 
between quarks are usually modeled as V(r) = A/r + Br, however it has been shown that 
A t dt = (A/r + Br)dt is not a solution of the Yang- Mills equations in the presence of a point 
charge [54]. The solutions (22), that we have found here, were applied to describe the energy 
spectra of quarkonia (charmonium and bottomonium) in [50], to predict the electric form 
factor, the magnetic moment and the root-mean-square radius of mesons in [55, 56, 57], and 
to study the chiral symmetry breaking in QCD [58]. 

It was shown that the confining solutions (22) satisfy the so-called Wilson confinement 
criterion [52]. This criterion is in essence the assertion that the so-called Wilson loop W(c) 
should be subjected to the area law for the confining gluonic field configuration. As a 
consequence, the latter law is equivalent to the fact that energy E(R) of the mentioned 
configuration (gluon condensate) is linearly increasing with R, a characteristic size of some 
volume V containing the condensate [52]. The evaluation of E(R) was carried out in [52] 
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a(r) 
a(r)0(r 



dr + B x , 

dr + Ai. 



(20) 
(21) 



using the Too - component of the energy-momentum tensor for a S'?7(3)-Yang-Mills field. 
However, formally E(R) diverges for everything R. Calculations carried out in [52] have 
considered the integral about the angle 9 inside the limits (9q,tt — 9 ), with the purpose of 
avoiding the divergence. Next we will consider a classical estimate of Wilson loop in the 
real-time formalism (r — » it) directly for the gluon condensate, i. e. 

W(R,T) =< e i9 $ A ^ >= e~ TV{R \ (24) 

and taking into account that the expectation value coincides with the evaluation of the 
integral in this approach. We obtain as result that 

//*2"7r 
A^dx^ = J {b x R + B 1 )TrAdip = 2vr(&i-R + BjTrA, (25) 

where Tr is the trace of the matrix. We observe from this result that the confinement 
potential obtained has the form V(R) = <jR + C. We note that this solution has the form of 
the confining Cornell potential [59] which has been used to describe experimental features 
of QCD. Lattice QCD simulations carried out have found the same kind of potential from 
Wilson loop [60, 61, 62]. 

C.2. Case of anti-de Sitter metric 

As a particular application of the (3 + 1) case solution, we will consider the anti-de Sitter 
metric given by 

ds 2 = (1 - Ar 2 /3)dt 2 - (1 - kr 2 /?>Y l dr 2 - r 2 (d9 2 + sin 2 9d<p 2 ). (26) 

As a(r) is the inverse of /3(r), then the Coulomb solution has not deformations respect to 
the flat space-time case, then the solution of the function f(r) is given by f(r) = ai/r + A±. 
On the other hand, the linear solution g(r) changes respect to the flat case. We get the 
solution explicitly in the two following situations A > and A < 0, so 

I h tanh- 1 + B u if A > 0, 

S(r) "\& lt an-(^) +Bl , i£A<0, (2?) 

The function g(r), for the limit cases |A| << 1 and r << 1, has the form g(r) ~ b±r + Bi, 
recovering the flat space-time case behavior. 

C.3. Case of Schwarzschild metric 

Another application for the (3 + 1) curved space-time solution is the Schwarzschild metric 
given by 

ds 2 = (1 - 2M/r)dt 2 - (1 - 2M/rY l dr 2 - r 2 (d9 2 + sm 2 9d<p 2 ). (28) 
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As in the last application, the Coulomb solution has not deformations respect to the flat 
space-time case, but the linear solution has. The function g(r) for this case is 

g (r) = b 1 (r + 2M\n\r-2M\) + B v (29) 

For the limit r >> M, this solution has the form g(r) ~ b±r + B±, recovering the flat 
space-time case behavior. 

We observe that the solutions (27) and (29) are confining solutions. From the analysis 
of the spectra of quarkonia in a spherical symmetry background, we can observe that the 
linear solution changes as a consequence of the curvature. This fact could generate a change 
in the mass spectrum of mesons. This change additionally could be generated by the the 
effects of the curvature in the Dirac equation. By this reason, the solutions (20) and (21) 
can represent a first step in the study of hadronic spectrum in a space-time curved. 

3.2. Case n = 2 

We consider now a curved space-time in (2 + 1) dimensions defined by the metrics 

ds 2 = g^dafdx" = a 2 {r)dt 2 - /3 2 (r)dr 2 - r 2 d8 2 . (30) 

For this case we assume that the connection (4) has the following form 

A = A M (r)dx" = A t (r)dt + A r (r)dr + A (r, 9)d9. (31) 

Using the gauge condition (8), we can obtain the following equation 

In the last equation, the solutions of the form Ag(r,d) = H(r)9 + G(r) can be discarded 
because they must be periodic in 9, i. e. Ag does not depend on 9. Thus, the equation (32) 
can be written as 

d r ( TO( 2 A ; (r) > | = 0. (33) 



The solution of this equation has the form A r = C^^. In these coordinates, the exterior 



/3(r) 

aLlWll 11CXD MIC 1VJ1111 T~L r — 

differential is written as 

d = O t dt + d r dr + d e d9. (34) 

If now we set C = and we substitute the functions A t = f(r)T and Aq = g(r)A in (31), 
where T and A are linear combinations of the group generators, and additionally we consider 
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the following relations 

r 

*(dtAdr) = — d9, 

a{r)B{r) 

, , ,„s a(r)8(r) , 
*(dt A d6) = v JHK ' dr, 
r 

OL\T ) 

*(dr A d0) = — ^rrft, 
rp(r) 

then the curvature is given by 

F = dA + #A A A = -d r f(r)Tdt A dr + d r g{r)Adr A dip + gf(r)g(r)[T, A]dt A dO, (35) 
and the Hodge star operator applied over the curvature can be written as 

„ rd r f(r)T ,„ a(r) d r q(r) . , a:(r)B(r)f(r)q(r) r ^ . . , , , 

*F= ,\\,\ d6+ „,\ yKJ Adt + g [T,A]dr. 36 

a(r)p(r) rpyr) r r 

The Yang-Mills equation for this case can be obtained from the expression (6) and has the 
form 

dr ( r d ;(^) A dr Add — d r ( a( - r ^r) A dt A dr = s ^?^± Tdr A dQ 
V a(r)3(r) J \ rB(r) J w a(r) 

of f(r)a(r) 2 a(r)B(r) rr ^ . . . . , ,„ f(r) 2 q(r)a(r)B(r) rr ^ , \ 

+ ^ 2 m m ; — i r, A , A cZr Ad6 — lA-UlAJ v r, A , r A . 

\ r r / 

(37) 

If the abelian condition is satisfied in the last equation, i. e. the commutator [r, A] = 0, 
then the solutions will be strongly restricted by this condition. Thus, we can obtain the two 
following independent equations 

dr(^d r g(r)) = 0, (38) 



rB(r) 

(T^TT^/tr)) = gStf^T, (39) 
\a{r)B(r) J a(0) 



which solutions are given by 



I r I^l dr + d2y (40) 

l^r^r) dr + di (41) 
9 r 



For a = (3 = 1, these solutions can be written as 

g(r) = d 2 + k 2 r 2 , (42) 
f(r) = di + kx logr. (43) 

We observe that the function f(r) has the form of the well known confining solution for 
the two dimensional problem [63]. Similarly, the function g(r) has a form of a confining 
function. It is possible to eliminate the constant d\ by means of a gauge transformation, so 
the potential is given by 

A^dx^ = hT log rdt + (d 2 + k 2 r 2 )Adtp, (44) 
and the 1-form vectorial field is 

A At = |A; 1 riogr,0,(^ + A; 2 r)A|. (45) 

In the last expression, the term kiTlogr is a potential equivalent to the Coulomb potential 
in (2 + 1) dimensions and additionally it corresponds to the potential which is obtained for 
an infinite line of charge in (3+1) dimensions. Additionally, the term k 2 r in the z axis is 
equivalent to a constant magnetic field of magnitude B\ = 2k 2 in (3 + 1) dimensions. We 
can see that if we perform a transformation to cartesian coordinates over this last term, we 
obtain its usual representation 

k 2 re v = -y r (- sin Wx + cos ipe y ) = -^-ye x + ^-xe y . (46) 

Finally, the term — can be interpreted as the potential produced by an infinite solenoid of 
radius R in (3 + 1) dimensions, with R — > 0. This last term can be written in its usual 
representation by means of a transformation to cartesian coordinates 

di A d 2 . ^ d 2 „ d 2 „ , > 

—e v = -(-sin pe x + cos (pe y ) = — -ye x + -^xe y . (47) 

This result is known as the Aharonov-Bohm (AB) potential, where the magnetic field within 
the solenoid is B 2 = d 2 . 

Let's consider now the integral that appears in the Wilson loop 

A^. (48) 

This integral can be performed in a circular path of radius R which is centered in the origin 
of the coordinates frame. Then we obtain that 

A^dx^ = (d 2 + k 2 R 2 )TrAdip = 2n(d 2 + k 2 R 2 )TrA. (49) 
Jo 1Q 



This result implies that for this approach the confinement potential in low dimensionality is 
given by V(R) = aR 2 + B. 

The result obtained for this case correspond to a U(l) abelian solution which has a 
confining behavior and can be considered as an application of the planar electrodynamics, 
i.e. electrodynamics in two spatial dimensions (QED 2 +i). QED 2+1 has been worked in many 
emblem quantum systems such as the quantum Hall effect, the theory of anyons and the 
relativistic quantum Hall effect [64, 65, 66]. During the last years the physics of graphene 
has attracted considerable interest, both theoretically and experimentally [67, 68]. Due to 
the low energy excitations of graphene can be described by a massless Dirac equation in two 
spacial dimensions, the curved graphene has been modeled by coupling the Dirac equation 
to the corresponding curved space [69]. In the approach presented in [70], gauge fields has 
been considered as external fields into the Dirac equation and it was possible to model some 
topological defects in the graphene. In connection with the results that we present in this 
work, the solutions that we obtain for the case n = 2 can be interpreted as a gauge field 
which is affected by the curvature of graphene. 

3.3. Case n = 1 

For this case we start from the metric given by 



where A a are the group generators. Each of these generators has associated a function f a (x). 
The gauge condition (8) implies that 



and then A x = C(3(x)/a(x). On this coordinates, the exterior differential is written as 



ds 2 = g iliV dx p 'dx v = a 2 {x)dt 2 — (3 2 (x)dx 2 . 



(50) 



We suppose that the connection has the following functional dependence 



A = A^(x)dx^ = A t (x)dt + A x (x)dx = X a f a (x)dt + A x (x)dx, 



(51) 




(52) 



d = d t dt + d x dx. 



(53) 



Fixing C = and using the fact that *(dt A dx) 



, it is possible to write 



Xgd x f a (x) 

a{x)f3{x) 



(54) 
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For this case we have that *(dt) = ^^dx and the expression (6) leads us to 

d x ( ^4M ) dx = g (^P-J b (x)[X a , X b }) dt + g5(x)^f\q a X a dx. (55) 
\ a{x)p(x) J \a{x)p{x) J a(x) 

A non-trivial solution from this equation can be obtained if the condition [A a , A&] = is 
satisfied, or alternatively X a d x f a (x) = d x f(x)T, where T represents a linear combination of 
the group generators. Thus, the last equation can be written as 

, f Xgd x f a (x) \ 0(0) 

d *{-^^)=- 95{x) W) qXa - (56) 

We solve this equation and obtain 

f a (x) = k a [ a{x)f3(x)dx + d a , (57) 



where k a and d a are constants. For a = (3 = 1 we can write f(x) = d + k\x\, i. e. we 
have obtained the solution of the flat space-time case. This last solution is invariant under 
a parity transformation and has the form of a linear confining solution. The form of this 
solution corresponds to the potential that results due to an infinite sheet of charge in (3 + 1) 
dimensions. 



4. Conclusions 

We have obtained some exact static solutions for the SU(N) Yang-Mills equations in a 
curved space-time of (n + 1) dimensions for the cases n = 1,2,3. For the (1 + 1) case, we 
have found that the solution for the temporary part can be written as A t = f(r)T. This 
solution is the most general static solution in a space-time in presence of a point charge. 
To be able to find analytic solutions in the (2 + 1) case, it was necessary to demand that 
the abelian condition given by [A, T] = were satisfied into the Yang-Mills equation for this 
case. For the cases (1 + 1) and (3 + 1), the abelian condition was naturally satisfied. We have 
presented in detail the solution for (3 + 1) curved space-time case and we have applied this 
solution to the anti-de Sitter and Schwarzschild cases. In both cases, the Coulomb solutions 
have not deformations respect to the flat space-time case, while in the linear solution exists 
deformation. We have assumed that these solutions can be considered as a first step in the 
study of the corrections on the spectra of quarkonia in a curved background. Although the 
energy diverges for the solutions that we have found, this fact is not a problem since this 
energy behavior was already known and it is not an impediment to make physical predictions 
in agreement with experimental data [50, 54, 55, 56, 57, 58]. Since the solutions that we have 
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found here are valid also for the group U(l), the case n = 2 is a description of the (2 + 1) 
electrodynamics in presence of a point charge. Finally, we have found that the solution for the 
case n — 1 is invariant under a parity transformation and has the form of a linear confining 
solution. As a perspective of this work, it would be interesting to understand the role of the 
confining solutions in a model of relativistic quark confinement in low dimensionality and 
the spectra of quarkonia in a curved background. 
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